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Abstract 


The n-Queens problem is to place n chess queens on an n by n chessboard so that no 
two queens are on the same row, column or diagonal. The n-Queens Completion problem is 
a variant, dating to 1850, in which some queens are already placed and the solver is asked 
to place the rest, if possible. We show that n-Queens Completion is both NP-Complete 
and #P-Complete. A corollary is that any non-attacking arrangement of queens can be 
included as a part of a solution to a larger n-Queens problem. We introduce generators of 
random instances for n-Queens Completion and the closely related Blocked n-Queens and 
Excluded Diagonals Problem. We describe three solvers for these problems, and empirically 
analyse the hardness of randomly generated instances. For Blocked n-Queens and the 
Excluded Diagonals Problem, we show the existence of a phase transition associated with 
hard instances as has been seen in other NP-Complete problems, but a natural generator for 
n-Queens Completion did not generate consistently hard instances. The significance of this 
work is that the n-Queens problem has been very widely used as a benchmark in Artificial 
Intelligence, but conclusions on it are often disputable because of the simple complexity of 
the decision problem. Our results give alternative benchmarks which are hard theoretically 
and empirically, but for which solving techniques designed for n-Queens need minimal or 
no change. 


1. Introduction 


The n-Queens problem is to place n chess queens on an n by n chessboard so that no 
two queens are on the same row, column or diagonal. This puzzle dates to 1848, and 
only two years later a variant was introduced by Nauck (1850) in which some number of 
queens are pre-placed and the solver is asked to place the rest, if possible. This is the n- 
Queens Completion problem and Figure 1 shows the first known instance studied. We will 
show that the n-Queens Completion problem is NP-Complete and #P-Complete, discuss 
solvers for the problem, and empirically analyse randomly generated instances. The n- 
Queens Completion problem may be one of the simplest NP-Complete problems to explain 
to people who understand the rules of chess. The problem is “Given an n x n chessboard 
on which some queens are already placed, can you place a queen in every remaining row so 
that no two queens attack each other?” 

The n-Queens problem has an extraordinary history for such an apparently unassum- 
ing problem, both generally and inside Artificial Intelligence. Formerly, and incorrectly, 
attributed to Gauss, the problem’s history was clarified by Campbell (1977). The 8-Queens 
problem was introduced by Bezzel (1848) and by Nauck (1850) (possibly independently). 
The latter publication attracted the interest of Gauss, who even made a small mistake in 
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Figure 1: This is the first published instance of the n-Queens Completion problem, by Nauck 
(1850). The reader may enjoy attempting to place 6 more queens on the chess- 
board so that no two queens attack each other. Is it possible? If so, how many 
different ways are there to do it? The answers to these questions are given below 
in Figure 12. 


studying the problem.! The generalisation to n-Queens has attracted the interest of many 
other mathematicians, whose results are surveyed by Bell and Stevens (2009). The first paper 
to describe backtracking search on computer, presented in 1958, does so for the n-Queens 
problem (Walker, 1960). Since then it has been used as an example and/or benchmark 
problem in many classic AI papers, for example at least six with more than 400 citations 
on Google Scholar at time of writing (Golomb & Baumert, 1965; Bitner & Reingold, 1975; 
Mackworth & Freuder, 1985; Minton, Johnston, Philips, & Laird, 1992; Selman, Levesque, 
& Mitchell, 1992; Crawford, Ginsberg, Luks, & Roy, 1996). 


The complexity of the n-Queens problem is often misunderstood. The decision problem 
is solvable in constant time since there is a solution for alln > 3 so is only NP-hard if P=NP. 
A witnessing solution can be constructed easily (Bell & Stevens, 2009) but note that the 
witness (a set of n queens) requires nlogn bits to specify but this is not polynomial in the 
size of the input, which is only logn bits. The n-Queens problem has often been шсоттес у 
called NP-hard, even in well-cited papers (Mandziuk, 1995; Martinjak & Golub, 2007; Shah- 
Hosseini, 2009; Nakaguchi, Kenya, & Tanaka, 1999, each with at least 29 citations). The 
counting version of the problem, i.e. to determine how many solutions to n-Queens there 
are, is sequence A000170 of the Online Encyclopedia of Integer Sequences (Sloane, 2016). 
The sequence is currently known only to n = 27, for which the number of solutions is more 
than 2.34 x 1017. No approach better than optimised exhaustive search seems to be known: 
e.g. the п = 27 total was counted using a massively parallel search using FPGAs (Preufer, 
2016). Hsiang, Hsu, and Shieh (2004) show that solving the n-Queens counting problem is 
“beyond the #P-class”. Bell and Stevens (2009) states that this means that there is no closed 
form expression in n for the number of solutions, but Chaiken, Hanusa, and Zaslavsky (2015) 





1. He reported finding 76 solutions but later realised that four of those were erroneous so had only 72. 
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claim to give one.? Cadoli and Schaerf (2006) studied the n-Queens Completion problem 
without studying its decision or counting complexity. The most closely related work to ours 
is by Martin (2007), who proves a rather different generalisation of n-Queens to be NP- 
complete, the key difference being that some squares can be marked as stopping attacks. 
In particular, this means that solutions to Martin’s problem need not be solutions to the 
n-Queens problem, unlike n-Queens Completion. 


We are contributing to a rich literature on the complexity of puzzles and games. One of 
the earliest results in the area is that generalised chess is EXPTIME-complete (Fraenkel & 
Lichtenstein, 1981). Amongst other games to have been proved NP-complete or harder are 
the card solitaire games Klondike (Longpré & McKenzie, 2009), Freecell (Helmert, 2003) and 
Black Hole (Gent, Jefferson, Kelsey, Lynce, Miguel, Nightingale, Smith, & Tarim, 2007), the 
Sudoku puzzle (Takayuki & Takahiro, 2003), video games like Pac-Man (Viglietta, 2014), 
and casual games such as Minesweeper (Kaye, 2000), Candy Crush Saga and Bejeweled 
(Walsh, 2014; Guala, Leucci, & Natale, 2014). Many other games are surveyed by Kendall, 
Parkes, and Spoerer (2008) and by Demaine and Hearn (2009). In each case, the complexity- 
theorist must define a generalised class of instances. Unlike some of the examples above, 
this step is completely natural for n-Queens Completion. 


Because of the ease of finding a solution, the n-Queens problem has been subject to re- 
peated controversy in AI over whether it should be used as a benchmark at all. For example 
a sequence of papers argued the point in the pages of SIGART Bulletin in the early 1990s 
(Sosic & Gu, 1990; Johnson, 1991; Bernhardsson, 1991; Gu, 1991; Valtorta, 1991), and then 
in 2014 the issue was raised again in a blog post by Smet (2014). We resolve this issue in 
the sense that, as an NP- and #P-Complete problem, n-Queens Completion does provide 
a valid benchmark problem. Similarly, the Quasigroup Completion Problem (which is to 
complete a partially filled latin square) is NP-Complete (Colbourn, 1984) and is a challeng- 
ing and popular benchmark (Gomes & Selman, 1997), whereas constructing a latin square 
from scratch is trivial. A very closely related problem, “Blocked n-Queens”, has previously 
been used for benchmarking without complexity guarantees (Namasivayam & Truszczynski, 
2009). Our results show as a corollary that Blocked n-Queens is NP-Complete and #P- 
Complete. We explore the practical difficulty of these problems, and a new variant, the 
Excluded Diagonals Problem. For Blocked n-Queens and the Excluded Diagonals Problem, 
we show the existence of a phase transition associated with hard instances, but we were not 
able to generate consistently hard instances for the n-Queens Completion problem. 


The remainder of this paper is structured as follows. We introduce a number of basic 
terms and concepts in Section 2. We give definitions of a sequence of problems in Sec- 
tion 3, where the first problem is n-Queens Completion, and the last is a variant of Boolean 
Satisfiability (SAT) that is NP- and #P-Complete, with the Excluded Diagonals Problem 
at an intermediate stage. The proof proceeds by a sequence of polynomial reductions in 
Section 4, starting with the last problem and ending with n-Queens Completion. Finally, 
we present an empirical study of random n-Queens Completion, Blocked n-Queens, and 
Excluded Diagonals Problems in Section 5. 





2. It is unclear to us if this is a mathematical dispute or simply a dispute on what it means to be a closed 
form expression, but in any case Chaiken et al.’s formula has not been used to extend knowledge of the 
number of solutions beyond n = 27. 
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2. Background and Definitions 


We will represent a queen as an ordered pair of integers. Since these represent coordinate 
positions in a grid, we will use vector notation where appropriate, although we will not use 
any significant vector algebra. Given a queen represented by an ordered pair (a, 8), the 
value а represents the queen’s column, and В its row on the chessboard. The values a + В 
and a — В represent the two diagonals the queen is on. It may not be intuitively obvious 
that chessboard diagonals correspond to sums and differences, but consider moving one 
square along the two orthogonal diagonals: in one direction the sum of the coordinates does 
not change, while in the other direction the difference does not change. To avoid possible 
confusion about which direction is which, we follow Bell and Stevens (2009) in simply calling 
a+ В the ‘sum-diagonal’ and a — В the ‘difference-diagonal’. 


Definition 1. A queen q is a pair of integers (a, В) with 0 < a, 0 < В. For the rest of 
these definitions we assume а = (a, В) and where appropriate a second queen q, = (a1, В1). 


® The “column of q”, written col(q), is a. The “row of а”, written row(q), is В. The 
“sum-diagonal of q” is a+ В or equivalently col(q) + row(q). The “difference-diagonal 
of q’ is a — В or equivalently col(q) — row(q). 


e The “size of а’, written ||q|| is defined by |9] = max(a, В) = max(col(q), row(q)). 
We say that а “fits on board size п” iff ||q|| < п. We extend this to sets of queens, so 


that ||Q]| = maxgeg((|q||). As normal, we write |Q| for the number of elements т Q. 


e We say that “q attacks q,” if they are not equal and any of the row/column/sum- 
diagonal/difference-diagonal is the same in а and ат; i.e. а attacks q, а а! and 
any ора = ол, P=fP},a+8=a,+)), ота-— В =а\ - В1. 


® Given an integer i, we write iq = (ia,iß) for the product of i and q. We write 


а + а! = (@+ «1,8 + 81). We extend these notations to sets of queens Q and Qı in 


. def г. def def 
the natural ways: iQ = {iq |q € О}, Q+qı = {a+a,|a€ Q}, О+О, = {q4+q, | 
Qe О, ат € Qi}. 


Problem 1 (n-Queens). PROBLEM: (n) where n is an integer > 1. SOLUTION: A set 
СО of queens which all fit on board size n such that: |Q| = п; and for any two distinct 
queens Q1,4» Е О, q, does not attack 42. If, additionally, for any two distinct queens 
(a1, 61), (a2, В2) Е О, both ay + 8, = a2+ By mod п and a, — В1 F a2 — Ba mod п, then Q 


is said to be а “modular” solution to the n-Queens problem. 


The modular n-Queens problem can be thought of as the n-Queens with the top and 
bottom edges of the board identified, and similarly the left and right edges identified, so the 
board forms a torus. Bell and Stevens (2009, Thm. 6) give the following result of Pölya. In 
terms of chessboards it says that we can obtain a solution to the n x m-Queens by taking a 
standard m-Queens solution and replacing each empty square with an n x n empty board, 
and replacing each queen with a copy of a solution to the modular n-Queens problem. It 
can be expressed in our terms as follows. 
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Theorem 2 (Pölya). If we have a modular solution Qı to the n-Queens problem, and any 
solution Qa to the m-Queens problem, then Qı + пФ is a solution to the п x m-Queens 
problem. 


In Theorem 2 the same modular solution is n-embedded at each one of m positions. We 
will need to embed distinct sets of queens to create one board so we define embedding at a 
single position. 


Definition 3. Given any queen q,, a set of queens Q2, and a natural number n (where 
|Qa|| < п), the n-embedding of Q2 at position qı is defined as nq, + Qo. 


3. A Sequence of Problems 


In order to prove that n-Queens Completion is both NP-Complete and #P-Complete, we will 
use a sequence of parsimonious reductions as described by Papadimitriou (1994, Chapter 
18). A parsimonious reduction maps instances of one problem to another such that the 
number of solutions is exactly preserved and the reduction can be computed in polynomial 
time. 

In this section we define n-Queens Completion and three increasingly general problems 
that generalise n-Queens Completion. We also define a restricted form of 1-in-3-SAT that 
is NP-Complete and #P-Complete. The proof (in the following section) is a sequence of 
four parsimonious reductions, starting with restricted 1-in-3-SAT and ending with n-Queens 
Completion. 

For each problem we define the parameters and the space of solutions. There is de- 
liberately no distinction made between the decision form of the problem (does a solution 
exist?) and the counting form (how many solutions are there?). When it is important, 
the distinction will usually be clear from context (for example, whether we are discussing 
NP-Completeness or #P-Completeness). 


Problem 2 (n-Queens Completion). PROBLEM: Ma = (n, P) where n is an integer and 
P is a set of queens that all fit on board size п, and no two queens in P have the same column 
or row as each other. SOLUTION: А set 52 of queens which is a solution to the n-Queens 
Problem for n and such that P C So. 


Note that this definition of n-Queens Completion does not require the input set P of 
queens to be mutually non-attacking, but does insist that they are all on separate rows and 
columns. This choice makes some cases in our later proofs easier, but makes no significant 
difference to solving complexity: if two queens in P attack each other on a diagonal, the 
answer ‘no’ can be given trivially. 

We next generalise n-Queens Completion by allowing restrictions that queens may only 
appear in a subset of the rows and columns (and must appear in each specified row and 
column). 


Problem 3. PROBLEM: M3 = (n, P,C,R) where (n, P) is an instance of Problem 2, and 
С, В are sets of integers such that |R| = |C| and for any queen (с, В) Е P, a € C and 
ВЯ В. SOLUTION: A set S3 of queens such that: |S3| = |C| + |P|; РС 53; no two queens 
in S3 attack each other; and for any pair (а, В) € 5з \ PwehweaeC, ВЕ В. 
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The reason that we exclude elements of P from С and В is to simplify later proofs. 
We next generalise Problem 3 by giving sets of sum-and difference-diagonals that must not 
appear. Unlike the row and column restrictions of Problem 3, there may be diagonals which 
are permitted but contain no queen in a given solution. 


Problem 4 (Excluded Diagonals). PROBLEM: Ма = (n,C,R,D~,D*), where M3 = 
(п, {}, С, В) is an instance of Problem 3 and D~, D* are sets of integers with D~ C {-(n— 
1),...,n-1}, Dt C {0,...,2n — 2}. SOLUTION: A set of queens 54 which is a solution to 
M3 and additionally: for any queen (a, В) Е Sı ше havea—8¢D-,a+8¢ Dt? 


Each of the problems above has a board size parameter n. We assume that the size of 
the encoding of a problem instance to a Turing machine input has size bounded below by n 
and polynomial in n. 

Our next problem class contains a set of instances of Problem 4, and ignores constraints 
between problems except for the sum-diagonals. Two queens placed within different sub- 
problems in the set are allowed to have the same row, column, or difference-diagonal, but 
not the same sum-diagonal. Thus we can compose instances of Problem 4 with control over 
their interactions. 


Problem 5. PROBLEM: A set Ms = {Мал | 0 <а< |M;|}, where each Мла is an instance 
of Problem 4 with D* = {}. SOLUTION: A set 55 = {S4a|0<a<|Ms|}, where each 54а 
is a solution to Мда, and additionally: for any {54а, San} С 55, and any (Qa, Bq) Е 54а, 
(ap, Вь) € Sap, we have aq + By F ay + By. 


We will use a restriction of the 1-in-3-SAT problem, in which all literals are positive, no 
variable occurs in the same clause twice, and no variable occurs in more than three clauses. 


Problem 6 (Restricted 1-in-3-SAT). PROBLEM: A pair Mg = (V,C) where C is a set (of 
clauses) such that each с € С is a set of size 3 of variables, с = {v;,vj;,uR}; and where 
У = {uv | de Е C.vu © с} is the set of all variables that are contained in any clause. 
Each variable v Е У occurs in at most three clauses in С. SOLUTION: A truth assignment 
56 : И — {true, false} such that for all c = {vj,v;, uz} Е С, 56 maps exactly one of vj, vj, UE 
to true and the other two to false. 





Theorem 4. Restricted 1-in-8-SAT is NP-Complete and #P-Complete. 


Proof. Porschen, Schmidt, Speckenmeyer, and Wotzlaw (2014, Lemma 4) prove that this 
problem is NP-Complete. The #P-Completeness of 1-in-3-SAT limited to positive literals 
(Monotone 1-in-3-SAT) is a corollary of Creignou and Hermann’s (1996) general results. This 
does not apply to the restriction on variable occurrences, so we must prove that Restricted 
1-in-3-SAT is #P-Complete. 

We will give a parsimonious reduction of Monotone 1-in-3-SAT to Restricted 1-in-3-SAT. 
Let nz, be the number of occurrences of variable x. We will create new variables that we 





3. Note that the sets D~,D* are diagonals that may not be used, while the sets В and C are rows and 
columns that must be used. This disparity is because we find ourselves specifying which columns are 
allowed (and implicitly the rest are disallowed), while usually specifying which diagonals are disallowed 
(and implicitly the rest are allowed). 
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force to be true and false, specifically {t?, ff, 5, 13.1, fg; |? =1-:. |2n„/3]}. For each i we 
add the three clauses 


{t, Ja) 1} {ti a1 f3it; {8 ЛЗ f3;-2} 


If some t? were false, then exactly three variable occurrences in the set of pairs { f3, 5, 11}, 
{f§,-1, f§} and {f§,, f},_2} would have to be true, which is impossible since each literal 
occurs twice. So these clauses force each t” to be true and each f” to be false, using only 
three occurrences of each t variable and two of each f. For each occurrence of variable x in 
the original clause set we introduce two new variables т; and x;. We create the clauses: 


{©1, 21, fT}; 122,82, fr}, | а Tr 


{is v2, erh {Z2, v3; fasah a) {n.-1 Ing» Е 


Since each f is false, the first line of clauses ensures that each x; and 2; take opposite values, 
while the second line of clauses ensures that all variables x; take the same value. Note that 
the clauses use at most two occurrences of each x; and 7;, leaving one x; to be used for a 
positive occurrence of x in the original clause set. We construct a clause set by replacing 
the it" (necessarily positive) occurrence of variable x in the original clause set with the 
variable x;, and then adding the clauses introduced above. A solution in which all x; are 
true corresponds with x being true in the original, while a solution in which all x; are false 
corresponds with x being false. This reduction is parsimonious because all values of x; and 
&; in the new problem are uniquely determined by the value of x in the original while all 
t variables are true and all f variables false. This uses at most three occurrences of each 
new variable and all are positive. Finally note that the translated problem has no variable 
repeated in a single clause: this is by inspection in all the new clauses, and for the translated 
versions of original clauses, note that every literal in every clause is replaced by a different 
variable. 














4. A Sequence of Reductions 


We will give a series of reductions from Problem 6 to Problem 2. Each reduction will be 
polynomial. Each will also be parsimonious, i.e. the number of solutions to the reduced prob- 
lem will be the same as the number of solutions to the original, thus simultaneously proving 
Problem 2 to be NP-complete and #P-complete. If we were interested in NP-completeness 
alone, we could use a simpler non-parsimonious reduction from 3-SAT to Problem 5 (as 
described in Section 4.4), but we are not aware of a way to simplify the sequence of reduc- 
tions from Problem 5 to Problem 2. Therefore the proof of NP-completeness alone would 
be almost as complex as the proof of both properties together. 


4.1 Filling in the Gaps: Reduction from Problem 3 to n-Queens Completion 


First we define sets of queens Qi9, An, By for use later in the reduction. 
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Figure 2: The solution Qi9 to the 19-Queens problem with an additional column and row 
marked in red. Two squares (17,19) and (19,8) on the additional column and row 
are marked with ‘m’, and are not on the same or adjacent diagonals to each other. 
Also emphasised are squares (0,0), (0,8) and (17,0). A number of squares are 
outlined to show that no queen is placed on the diagonal (0,0) to (18,18) or either 
adjacent diagonal and that no queen is on the diagonals or adjacent diagonals of 


squares (17,19) and (19,8). 


Definition 5. We define the following sets of queens, illustrated in Figures 2 and 3. 
Qio = 00,8), (1,12), (2,17), (3, 7), (4, 10), (5, 15), (6, 1), (7, 16), (8, 13), (9,2), 
(10, 6), (11, 18), (12, 3), (13, 5), (14, 11), (15, 9), (16, 14), (17, 0), (18, 4)} 
An = {(i,2imodn)|0<i<n} 
В, = {(2imodn,i)|0<i<n} 


Next we observe some properties of Qig9 that will be used in later proofs. 


Observation 6. Q19 is a solution to the 19-Queens problem. It has no queen on difference- 
diagonals 1, 0 от -1. For the queens (0,y),(x,0), consider the pairs (19, у) and (x,19): that 
is, in Qig consider (19,8) and (17,19) (marked with ‘m’ in Figure 2). Their sum-/difference- 
diagonals are not within 1 of each other, (0,0) or any queen in Qig. 


All these properties can be checked by inspection. They can also be visually seen in 
Figure 2. Qi9 was found by a search using Savile Row (Nightingale, Akgiin, Gent, Jefferson, 
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Figure 3: Two modular solutions А1з (left) and By3 (right) to the 13-Queens problem as 
defined in Definition 5 with each queen in А1з at (1,25 mod 13) and in Bj3 at 
(2i mod 13,72) (right.) The first seven columns on the left have no two queens 
in adjacent rows, while the first seven rows on the right have no two queens in 
adjacent columns. The two solutions are mirror images along the diagonal (0,0)- 
(12,12), and also use exactly the same sets of both sum- and difference-diagonals. 


Miguel, & Spracklen, 2017) and Minion (Gent, Jefferson, & Miguel, 2006) from a model 
available online.* Further Savile Row/Minion searches show that there is no n < 19 for 
which a solution to n-Queens can have all the properties of Observation 6. 


Lemma 7. Jfn =1 mod 6 orn=5 mod 6 then Ay and By, are modular solutions to the 
n-Queens problem, and use the same sets of sum-/difference-diagonals as each other. 


Proof. Bell and Stevens (2009, Thm. 9) prove that Ap is a modular solution to the n-queens 
problem for п = 1 mod6 or п = 5 mod6. As By, is symmetrically equivalent to Ay, 
(reflected about the longest difference-diagonal), this implies В» is modular as well. 
Because (x,y) € An iff (y, x) € By, the sets of sum-diagonals in A, and В», are the same. 
As n is odd, the difference-diagonals of A, are {i — 2i| i € {1.454}} U {i- (ii - o lie 
nt) п} }, which simplifies to {-i | ie {1.252} U!n-i lie {44 n}} = we SZ. 
Ber (x,y) Е An iff (y,x) Е В; the diferenos diagonals of A, are the negation of the 
difference-diagonals of Ви, therefore the difference-diagonals of A, and B, are the same. 




















Lemma 8 (Rivin and Zabih). Suppose we have a solution Q to the n-Queens problem, a 
subset of queens PC Q, and а set of queens P, which uses the same set of rows, columns, 
sum-diagonals and difference-diagonals as P. Then (Q \ P)U P, is also a solution to the 
n-Queens problem. 


Proof. This result is due to Rivin and Zabih (1992, Thm. 2.4). Since P; uses the same rows, 
columns and diagonals as P, no queen in © \ P can attack any queen in P,. 

















4. https://ipg.host.cs.st-andrews.ac.uk/n-queens-completion-experiments-code.tgz 
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The following lemma proves that if two queens are not in the same row, not in the same 
column, not on the same diagonal, and not on two adjacent diagonals in the same direction, 
then: if we multiply these two queens by n and add to each one any queen that fits on board 
size n then the resulting two queens will not attack each other. 


Lemma 9. Consider any two queens q,,Q 2 which both fit on board size п, and two queens 
93,94 of any size. Then: if the columns (respectively rows) of аз, ал are different, then 
the columns (rows) of а: + паз and 4> + пад are different. Also: if the sum-diagonals 
(respectively difference-diagonals) of аз, ал differ by more than one, then the sum-diagonals 
(difference-diagonals) of чт + паз and 42 + пад are different. If all the just-mentioned 
properties hold for 4з, 94 then q, + паз does not attack qo + пад. 


Proof. Let each а; = (a;,ß;). For columns, suppose аз # a4. Since 0 < aj,a2q < п, 
ол + паз # а2 + пад, so the columns of q; + паз and 42 + пад are different. The case of 
rows is similar. For sum-diagonals, suppose that |(a3 + 83) — (a4 + В4)| > 2. The range 
of sum-diagonals on a board of size п is from 0 to 2n — 2, зо |(aı + 61) — (a2 + Ba)| < 2n. 
But |n(a3 + Вз) - n(a4+ 64)| > 2п. Hence the sum-diagonals of q, + паз and 42 + nq4 are 
different. The case of difference-diagonals is similar, since their range is —(n — 1) ton—1, 
so |(aı — Ba) — (a1 — В5)| < т. If all these cases apply then а; + паз does not attack 
Qo + пад, by definition of not attacking. 














We can now provide a reduction from Problem 3 to Problem 2. This reduction takes 
an instance of Problem 3 and n’-embeds it at location (0,0) of a large board based on 
Figure 2 (where n’ is defined as part of the construction). The other queens in Figure 2 
are also replaced with chessboards of queens, which ensure the problem is a valid instance 
of Problem 2 and block the rows and columns required of the copy of Problem 3 in the 
bottom-left corner. To solve the constructed instance of Problem 2, queens must be placed 
within the embedded instance of Problem 3 and nowhere else, hence a solution of Problem 2 
maps back to Problem 3 very easily. 


Definition 10 (Reduction from Problem 3 to Problem 2). Given an instance of Problem 3, 
M3 = (ns, В, С, В), we construct an instance of Problem 2, Ma = (na, P2) as follows: 


e If a pair of queens {41,9} С P3 attack each other, then we simply drop В and C and 
set ng = пз, Pa = P3. Otherwise: 


1. If ng mod 3 = 2 then set п’ = 2n3 + 1, otherwise set п’ = 2n3 — 1, so Lemma 7 may 
be applied. 


2. Form the set of queens Qı = Ay +n’(Qı9 \{(17,0)}). Set Q2 = Qi U (By + n'(17, 0)). 
That is, n’-embed a copy of A, at position а for every q € О19, except for queen 
(17,0), for which we n’-embed a copy of Byr. 


3. Set С’ = СЦ {col(q) |q € В}, В' = RU {row(q) |qe P3}. Set Q3 = {а | qe 
05, со(а) ¢ C’,row(q) ¢ В}. That is, we remove queens from @5 if they are in one 
of the columns in С, rows in А, or one of the same rows or columns аз a queen in P3. 
This can only remove queens from the copy of A, that is n’-embedded at (0, 8) or the 
copy of B, that is n’-embedded at (17,0) in Qo. 
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4. Put the columns т С’ т ascending order ag <aı <... < Qic—1- Form the set of 
queens Qo = {(i,2a;) | 0 < i < |С'|}. Similarly order the rows By < By <... < 
Bir in Ri and set Qr = {(28;,%) |0 <4 < [В]. 


5. Set no = 19n’ + |C’|, Po = Q3 U (Qc + n’(19,8)) U (Qr + n'(17, 19)) U P3. 


We will now prove that the construction in Definition 10 is correct, and explicitly show 
how to map from solutions of Problem 2 to Problem 3 and vice-versa. 


Lemma 11. If 52 is a solution to Ma, then S3 = {а | аЕ 52 & ||q|| < пз} is a solution to 
M3. 


Proof. No queens in $3 attack each other since 52 is a solution. Certainly Рз С Pa by 
construction and P, С 52 since 55 is a solution. Since ||P3|| < пз this gives us P3 С 53. 
It remains to show that |53| = |C| + |P3|, and for every а Е $3 \ P3, both col(q) € С and 
row(q) € В. For a given column a € С, we have a < пз, and there must Бе qq € 52 with 
col(qq) = a, and this has row(qq) = В + jn’ for some 0 < В <n’ and 0 < 7 19. We 
will show that 7 = 0 and 6 € R. The only possible nonzero value of 7 is 8, since all other 
rows are excluded by copies of Ayr: the possible values of В for 7 = 8 are the rows occupied 
by queens deleted at step 3, ie. the rows {2a | a € C’} by construction of А». But each 
of these is excluded by one of the queens in Ос + n’(19,8) added at step 5, specifically by 
(i, 2a;) + n’(19,8) if a; is the i” element in С’ in ascending order. So we have j = 0. To 
check В € В, first note that all rows outside В’ are excluded by queens from В», + n’(17,0), 
and rows from queens in P3 are excluded by P3 С Pa. Зо ВЕ R’ \ {row(q) | q € P3} = В. 
That is, we have shown that for all columns a € С there is a queen (a,3) € 53 with 
ВЕ R. This gives |.$3| > |C| + |P3|. But n’ — (|C| + |P3|) possible columns are ruled out 
by queens remaining from Ay + n’(0,8), so |S3| < n’ — (n’ — (|C| + | P3|)) = |C| + |. So 
|S3| = |C| + |P3| as required, and no queens in $3 \ P3 can use a column not in С’ or row 
not in R. 














Lemma 12. If 53 is a solution to M3, then 52 = 53 Ч Pa is a solution to Ma. 


Proof. We note that 55 = 530 (Qc + п (19, 8)) U (Qr+n’(17,19)) U Оз. We will prove that 
no two queens within one of these four subsets attack each other, and that no queen in one 
subset attacks any queen in the other subsets. 

First we consider attacks within one of the subsets. $3 is а solution to M3, so no two 
queens in it attack each other. By Lemma 7, A, is a modular solution to n’-Queens, so by 
Theorem 2, Qı is a subset of a solution to 19n’-Queens (Qı U Ay + n’(17,0) is a solution 
to 19n’-Queens). Lemma 7 applies since n’ is constructed to be = 1 or = 5 mod 6, and 
with Lemma 8 shows that Q2 is a solution to 19n’-Queens, and Оз С О5, so no queens in 
Оз can attack each other. Attacks within Ос + n’(19,8) are equivalent to attacks within 
Qc. The queens in Qe are (i,2a;), with 1 = 0,1,2,... and a; strictly increasing with i. 
Queens within Qc are trivially on different rows and columns; the sum-diagonal 7 + 2a; 
is strictly increasing with 7; and finally 22541 — 2a; > 2 so i+ 1 — 2aj41 < i — 2a;, so 
the difference-diagonals are strictly decreasing with 7. So no queens within Qc attack each 
other. The case of Qr + п (17, 19) is similar, with the queens being {(28,,1) |0 <i < |R’|} 
and В; strictly increasing. 
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Now we consider attacks between queens in different subsets. By the properties of Qi9 
in Observation 6, and from Lemma 9, no queen in Ос + п’ (19,8) attacks any queen in the 
rest of 52 except possibly queens in Ay + n’(0,8). But any queen in A, + n’(0,8) on the 
same row as a queen in Ос + n’(19,8) is not in Оз by step 3. The case of Ов + n’(17, 19) 
is similar. The only remaining case is a queen in 53 attacking a queen in Q3. Again, we 
use Observation 6 and Lemma 9 to show that no queen in Sy can be on the same diagonal 
as any queen in Оз. Finally, all queens in $3 must be in a column in С’ and a row in В’, 
but the queens in those columns and rows were removed at step 3 to form Q3, so no attack 
along a row or column is possible. 














Theorem 13. The reduction of Definition 10 is a parsimonious reduction from Problem 3 
to Problem 2. 


Proof. From Lemma 11 and Lemma 12, and the fact that the maps defined in those lemmas 
are inverses of each other, since if 92 = $3 U P2, we have $3 = {q|q € S2 & ||q|| < n3}. 














4.1.1 INCLUDING ARBITRARY n-QUEENS PARTIAL SOLUTIONS 


Lemma 12 has simple corollaries that show any set of nonattacking queens can be found as 
a part of a solution to a larger n-Queens problem. While a brief detour from our complexity 
results, we find the result attractive and it may help illuminate the structure of the n-Queens 
problem. For example, it shows that no local arrangement of nonattacking queens can cause 
global unsolvability for all n. First we define inclusion of a set of queens P within a larger 
set Q. Inclusion is an embedding with the additional requirement that no queens that are 
not part of P occur in that subgrid of Q. 


Definition 14. If P is a set of queens with || P|| < п, we say that P is n-included at location 
(т, у) in a set of queens Q if P+(x,y) С О and there does not exist a queen q with ||q|| < п 
such that q ¢ P andq+(a,y) Е О. 


Corollary 15. For any set of queens P, where n = ||P||, with no two queens in P attacking 
each other, P may be n-included at (0,0) in a solution of the na-Queens problem for some 
na. 


Proof. Set M3 = (n, P,{},{}), and construct an n-Queens Completion problem Ma = 
(na, Pz) as in Definition 10. P is a solution of Ms, therefore (by Lemma 12) Ps is а so- 
lution of Ma that n-includes P at position (0,0). All solutions of Ma = (na, P2) are also 
solutions of the n-Queens problem Mı = (na) so the corollary holds. 














Not only can we include P at (0,0), it can also be included at an arbitrary position on a 
larger grid simply by applying Corollary 15 to Р+ (1,1) for any position (7,7). Furthermore, 
we can include as many copies of P as we want in a single solution to a larger n-Queens, 
and do so in exponentially many different ways. First, consider the four sets P + (0,0), 
P+(0,1), P+ (1,0), P+ (1,1). Assuming P is non-empty, these four sets are distinct and 
all include P. Second, consider the set of К queens K = {(3i,i) |0 <1< К}. The queens in 
K are all on different rows and columns and their sum-diagonals all differ by more than one 
from each other, as do their difference-diagonals. By Lemma 9 we can (n + 1)-include any 
one of the four sets that include P in positions in (п + 1)K, giving a larger set of queens 
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P'. There are 4" ways of constructing P’. By applying Corollary 15 to P’ we have 4* ways 
of including k copies of P in a solution to a larger n-Queens problem. 


4.2 Ruling out Diagonals: Reduction from Problem 4 to Problem 3 

Definition 16 (Reduction from Problem 4 to Problem 3). For the instance of Problem 4 
Ма = (n,C,R,D”,D*), set Ор- = {(6n — 3 -d,3n—1-2d) | de О}, set Qp+ = 
{(2n—-2-d,n+2d) |de Dt}, and set МЗ = (7n-3,Qp- UQp+, {c+3n—2| cE С}, В). 





The reduction is illustrated in Figure 4. The key point of this construction is that the 
difference-diagonal of (6n — 3 — d,3n — 1 — 2d) is 3n — 2 + 4, and the sum-diagonal of 
(2n — 2—d,n+ 2d) is 3n — 2+.d. The offset of 3n — 2 in both cases explains why we shift 
the set C by 3n — 2 in M3. 

First we confirm that M3 is a well-defined instance of Problem 3, assuming My is itself 
well-defined. The range of possible values of D~ is —(n—1)...n—1, and hence the extreme 
possible placements for queens in Qp- are (7n — 4,5n — 3) and (5n — 2, п +1); D* ranges 
from 0 to 2n — 2, so the extreme possible placements for queens in Qp+ are (0,5n — 4) and 
(2n —2,n). These certainly have no overlap with the sets Rand {с+3п-—2 | сЕ С}, and all 
fit on the board size 7n — 3 of M3. No queen in Qp- UQp+ is on the same row or column as 





any other: the only non-trivial observation necessary is that if n is odd (even), then queens 
in Qp- are all on even (odd) rows and Ор- on odd (even) rows. So M3 is well-defined. 


Lemma 17. If S4 is a solution to Ma, then S3 = (54 + (3n — 2,0)) UQp- U Qp+ ва 
solution to M3. 


Proof. Within each subset of $3, no pair of queens attack each other. In the case of 54 + 
(3n — 2,0) this is because 54 is a solution, and in the other subsets by construction. We 
also noted above that queens in Qp- and Qp+ are on different rows from each other, and 
their columns and diagonals are far apart. So we only need to consider potential attacks 
between Qp- UQp+ and 54 + (3n — 2,0). But if (a, В) Е 54 then a— 6 ¢ D™ and so 
(6n — 3 — (a — B),3n — 1- 2(a - В)) ¢ Qp-, therefore (a + 3n — 2, В) is not attacked by 
any queen in Qp-. Similar reasoning applies for Qp+. Finally we note that constraints on 
columns and rows in 53 are obeyed because they are in МА. 














Lemma 18. If Sy is a solution to Ms, then 54 = {(a — Зп +2, В) | (a, В) € 53,8 < n} is 
a solution to Ma. 


Proof. No two queens in 54 attack each other since 553 is a solution. They all come from the 
first n rows of 53, so the range of a — 3n +2 is 0 to n—1, so ||S3|| < п. The constraints on 
membership of C and В are satisfied since they are in M3. For any de D7, a queen in Qp- 
is on difference-diagonal 3n — 2+ 4, so no queen in $3 is on that diagonal, hence no queen 
in 54 is on difference-diagonal 4. Similarly for diagonals in D*, so the diagonal constraints 
in Ma are obeyed. 














Theorem 19. The reduction of Definition 16 is a parsimonious reduction from Problem 4 
to Problem 8. 


Proof. From Lemma 17 and Lemma 18, and the fact that the maps defined in those lemmas 
are inverses of each other. 














827 


GENT, JEFFERSON, & NIGHTINGALE 





\ 
\ 
\ 
\ 





_ 





пен 
N 
ни, 
ma 


m 
> 


un ия “a 
a me : 
a Y а 


= 
w 





BRR 
[> 


OrRFN WHF UU DN O O 














Figure 4: Ruling out the diagonals of the 4 x 4 square in bottom centre of the grid by 
Definition 16. A пх п square has 2n — 1 diagonals in each direction. In this 
example we see 7 queens from (0,16) to (6,4) ruling out the sum-diagonals, and 
7 queens from (18,5) to (24,17) ruling out the difference-diagonals. For example, 
the queen at (4,8) rules out the sum-diagonal 12, while the queen at (21,11) rules 
out the difference-diagonal 10. Note that none of these queens attack any squares 
in the 4 x 4 square except those in the intended diagonal, and also that no two 
of the 14 queens attack each other. Therefore we can use any subset of these 14 
queens to rule out exactly and only a selected set of diagonals in the 4 x 4 square. 


4.3 Arrangement of Gadgets: Reduction from Problem 5 to Problem 4 


The reduction from Problem 5 to Problem 4 is relatively straightforward. Since an instance 
of Problem 5 is made up of multiple instances of Problem 4 with interacting sum-diagonals, 
it is enough to lay those instances out on the same long diagonal of a larger grid, separated 
enough that their rows, columns, and difference-diagonals do not have any effect on each 
other. Apart from combining the separate subproblems correctly, we must ensure that no 
queens can be wrongly placed in any square off the intended long diagonal. To do this, we 
exclude all sum-diagonals which do not intersect with any of the instances. Figure 5 gives a 
high-level illustration of the reduction. 


Definition 20 (Reduction from Problem 5 to Problem 4). Let Мь = {Мда | а =0,1,2...k— 
1} be an instance of Problem 5 where К = |Ms| and each Maa = (Na, Ca, Ra, Dz, {}). 
first set п’ = max{ng |0<a<k}. We will place each gadget My at position (aa, Ва) in a 
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[] Instance of Problem 4 


[а Ruled out Бу row or 
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Ruled out by sum- 
diagonal constraint 











Figure 5: A schematic view of the reduction from Problem 5 to Problem 4 in Definition 20 
for an instance of Problem 5 containing 5 instances MO to М4 of Problem 4. 
Each square shown represents an n’ x n’ cell, where n’ is the size of the largest 
subproblem. The instances МО to M4 are placed on the long sum-diagonal, from 
n’(0,8) for МО to n’(8,0) for M4, indicated by red outlines. Green areas indicate 
that all cells on odd-numbered rows and odd-numbered columns are disallowed. 
Sum-diagonals that are not on the long diagonal are excluded, and this is indicated 
by the grey areas. The additional rows, columns and diagonals excluded within 
each red cell are not shown. 


nxn grid and define the instance Problem 4 as My = (n,C,R,D~,D*), where: 


Qa 


Ba 


n 


Dt 


2n'a 

2n'(k —1—a) 

т (2k — 1) 

U {c+a.|ce Ca} 
0<a<k 

U ee 
0<a<k 

U {d+ oq=—8,|¢eD,} 
0<a<k 


{0,1,2,...2n’(k —1) - 1} U {2n'k, 2n’k +1,...An’k — 2} 


The definition of D* is very different from the others, which are simple aggregations. 
The role of D* is to exclude all diagonals which are outside the range [?2n’(k-1),2n’k). The 
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reason for this is to exclude queens at positions which cannot be excluded by either C or В 
but which cannot be allowed: these are seen as grey cells in Figure 5. The well-definedness 
of M, is straightforward, given the well-definedness of each component Мл. 


Lemma 21. 155 = {54а | О < a < К} is a solution to Ms, then 54 = Uo<acr (aa + 
(ма, Ва) is a solution to Ma. 5 


Proof. The queen q, = (a, В) € 54а maps to Qu. = (а а, В- В.) € 54. First we show no 
two queens in 54 attack each other. Two queens in 54 mapped from the same S4,. both have 
the same offset (аа, 3,), so cannot attack each other since no queens in 54а do. For queens 
not from the same 54а, first consider columns, rows and difference-diagonals, and queens 
a Е 54а, Gy € Sa, а # b. Lemma 9 applies with (in terms of that lemma) n = п’, qı = 4.; 
G2 = 9; 93 = (2a,2(k —1—a)), 94 = (2b,2(k — 1 — b)), so q, and q can only possibly 
attack each other along sum-diagonals. Note that the sum-diagonal for q, is a+ В, while 
d4q Ваз sum-diagonal си + оз + В+ В. =at 6 +2n'a+ 2n'(k—-1—a) =a+6+2n'(k—1). 
Thus the sum-diagonal in $4 is the constant 2n’(k — 1) plus the value of the sum-diagonal 
in 54а. So two sum-diagonals in 54 are equal iff the corresponding sum-diagonals are equal 
in $5. Since 55 is a solution, that is not the case, so no two sum-diagonals in 54 are equal. 
Hence no two queens in 54 attack each other. 

We must also show that all restrictions based on C, В, D~, D* are obeyed by Qua. For 
С note that a Е Cg since 55 solves Ms, and зо а + аа € С, as required. The case of В is 
similar. The difference-diagonal of qy, is (a+ аа) — (8+ 8,). But a— ВЯ Dz because 
Dz represents the forbidden difference-diagonals of Ma... Therefore (a + aq) — (8+ 8,) ¢ 
{а+а. - 6, | de D7}. We need to check that we cannot have (a + aq) — (В+ В.) € 
{d+ ap — В, | d€ D, } for some b ¥ a, i.e. the difference (a + aq) — (8 + Bg) — (d+ an — By) 
cannot be zero. Rearranging we have a — В — d+2n/(a—-b-(k-1-a)+(k-1-b)) 
=a—8—d-+4n'(a—b). а- В -а| < 2n’ and |An’(a — b)| > An’, so the sum cannot 
be zero. Therefore the difference-diagonal of Чл и is not in О”. The sum-diagonal of qy „ 
is (@+ Qa) + (8+ By) = a+ та+ (B+ 2n'(kK-—1-—a)) =a+6+2n'(k—1). But 
0<a+ 6 < 2(п' — 1) so 2n’(k-1) < a+ В+ 2n'(k — 1) < 2n’k, so the sum-diagonal of 
iq is not in Dt. 

















Lemma 22. If Sı is a solution to Ma, then 55 = {54. | 0 <a <k} where 54а = {(a, В) | 
0<a,8 <n’, ((а, В) + (aa, 8q)) Е 54} is a solution to Ms. 


Proof. The key point is to show that all queens q € 54 are of the form q = (a+ аа, В + Ва) 
with 0 < a,6 <n’, for some a. The construction of C and В gives а = (а + аа, В + ß}) 
with 0 < a,8 < n’, but we have to show that a = 6. The sum-diagonal а of а is а = 
аа. +В+ 6, =at+ В+ 2n'a+2n'(k-1—6) = а В+ 2n'(k-1+(a—b)). But 
0 < а-+ В < 2n'— 2, зо На < b then а < 2n'(k — 1) - 1, while if > b then d > 2n’k. 
In either case 4 Е Dt, so we must have а = b. All of a € Cy, ВЕ В, а-- B ¢ р: 
now follow immediately from the corresponding facts for S4. We have ||S4,q|| < па because 
each max(C, U Ra) < na. Certainly no two queens in any 54а can attack each other since 
none in S4 do. Finally, no two queens in 54а, Sa, can be on the same sum-diagonal since 
the corresponding queens in $4 are not and the sum-diagonals in 54 are the constant value 
2n’(k — 1) greater than the sum-diagonals in each of 54а, 54. 
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Theorem 23. The reduction of Definition 20 is a parsimonious reduction from Problem 5 
to Problem 4. 


Proof. From Lemma 21 and Lemma 22, and the fact that the maps defined in those lemmas 














are inverses of each other. 


4.4 Gadgets for Clauses and Variables: Reduction from 1-in-3-SAT to 
Problem 5 


For an instance Mg = (V,C) of Restricted 1-in-3-SAT, we will construct an instance M5 = 
{Ga} of Problem 5. Recall that each Ga = (па, Ca, Ва, Dz, {}) is an instance of Problem 4.° 
We assume that У = {м | i = 0,1,...|V| —1} and С = {c; |i = 0,1,...|C| -1}. The 
ordering of both variables and clauses is arbitrary but must be fixed before the reduction is 
performed. We introduce the notation G„®6., which can be thought of as shifting Gg by да 
places along the x-axis. 

Definition 24. For an integer да > 0 and Ga = (Na, Ca, Ra, Dz, {}), we define Ga®da ae 
(na + да, {C+ ba | СЕ Ca}, Ra, {d+ 5a |de Dz}, {}). 


Lemma 25. A set of queens Оа is a solution to Ga iff Qa + (6a,0) is a solution to Ga®da 


Proof. The addition of (64,0) does not affect attacks between two queens, so we only have 
to consider the other properties of Qa + (да, 0), and each is easy. Оз fits on board size nq iff 
Qa + (6a, 0) fits on size na + ||(da,0)|| = па + да. For queen (a, В) € Qa, we have a Е C, Ш 
а-+ да = col((a@, В) + (да,0)) Е {e+4a | сЕ Ca}, and similarly for Ва. Finally, (a-ß) ¢ D7 
iff col((a, В) + (da, 0)) — row((a, В) + (4a,0)) =a—-—B+a ¢ {d+6.| de Dz}. 














С0: A gadget for а SAT variable. A copy of the same gadget G® is used for each SAT 
variable in a problem. Each copy of С° has two solutions, and we will arrange that one 
solution will indicate that the corresponding SAT variable should be true, while the other 
solution will indicate that the SAT variable should be false. The gadget G® is illustrated in 
Figure 6 and defined in Definition 26. G® has only two solutions, as can be checked by a very 
small search. The two solutions are {(0, 0), (1,3), (3, 2)}, which we call the ‘red solution’, and 
{(0, 2), (1,0), (3, 3)}, the ‘blue solution’. These use disjoint sets of sum-diagonals. These are 
respectively {0,4,5} (the ‘red diagonals’) and {1,2,6} (the ‘blue diagonals’). The gadget 
СО was found by a search using Savile Row (Nightingale et al., 2017) and Minion (Gent 
et al., 2006) from a model available online.® 

For a variable v; in a clause c we will use one red diagonal and one blue diagonal from 
the #” copy of G®. We will write а; . for the red diagonal being used, and dp for the blue 
diagonal, and d?, for a diagonal of indeterminate colour (ie. x € {b,r}). Note that in 
Restricted 1-in-3-SAT, v; occurs in at most three clauses, and no more than once in any 
clause. Because of this, for any given variable v, the value d,, . determines the values с and 
x. 





5. We use the letter G instead of Ma for gadgets in Ms to avoid proliferation of subscripts and to emphasise 
that these are gadgets, which will be used to encode different parts of a 1-in-3-SAT instance. 
6. https://ipg.host.cs.st-andrews.ac.uk/n-queens-completion-experiments-code.tgz 
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Figure 6: Illustration of the 4 x 4 gadget G°. The first diagram shows a board with arrows 
barring the second row and third column. The second and third boards show the 
only two ways that three queens can be placed in the remaining squares with no 
two queens attacking each other. The first, red, solution will be used to indicate 
that the corresponding SAT variable is true, with the second, blue, solution for 
false. The final board shows that the red and blue solutions occupy disjoint 
sets of sum-diagonals. Each square is labelled with the sum-diagonal it is on, with 
diagonals 0, 4, and 5 being red and 1, 2, and 6 being blue. The other sum-diagonal 
(3) is not used in either solution. 


Definition 26 (Gadget С°). We define С° = (4, {0, 1,3}, {0, 2,3}, {}, {}). For a variable 
v; where0 <i < |V|, we define dy, ae 214. (We sometimes write 6; for dy, to avoid the 
double subscript.) Each variable v; is associated with one copy of G° named G®, which is 
defined as follows: С = С0Фбь,. 


We allocate red and blue diagonals as follows. If v; occurs in clauses cı <c2 < с3, then 


d d d d 

for cı we set d;.. a 0, a? ., = 1; for co we set d;., = 4,4, = 2; and for сз ше set 
d d . В 

а; og = 5, = 6. If vu; occurs in fewer than three clauses we simply define the relevant 


subset of the dj... 


G!, (2, СЗ: Three gadgets for a 1-in-3-SAT clause. Now we present the gadgets 
associated with a clause с = {v;,vj, vg}, namely Gt, G?, СЗ. GI will ensure that at least 
one of the G®, Go and G® takes its red solution, while С? and G3 will each ensure that at 
least one of them (different in each case) takes its blue solution. Combined they will force 
exactly one of G9, Go and a to take its red solution: thus variables assigned true in SAT 


will correspond precisely with СО gadgets using their red solutions. The gadgets GL, G? 
and СЗ will vary from clause to clause, unlike G°, where the gadgets are the same for each 


variable. For a variable v € {v;,v;, vg} and colour x € {b,r} we define со = dy + 4, с: all 
columns in GL, G?, G3 will be of this form. The rows in G!, G2, G3 will also be of standard 
form, drawn from {7a | a € 0,1,2} = {0, 7, 14}. 


Definition 27 (Gadget G!). Suppose we have a clause с = {v;,0;, Un} of Restricted 1-in-3- 
SAT where 1 < 1 < К, so also 6; < 6; < дк. 


e Gi 2 (col? + 1, Си, Ri, D7, {}) where: Cy = {col?, col’, colk}, Rı = {0,7,14}, and 
Dy = {col? — 7, со — 14}. 
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Key 


[] Possible position of 
queen 


[] Ruled out by difference- 
diagonal constraint 





b b 
i ve cals, Ruled out by row or 
column constraint 


Figure 7: Illustration of Gl showing the seven positions where queens may be placed. 


Key 


Possible position of 
queen 


Ruled out by row or 
column constraint 


Interaction on sum- 
diagonal between two 
gadgets 





Figure 8: Illustration of С? and СЗ and their interaction through one sum-diagonal (col +0). 


Gadget G! is illustrated in Figure 7. There are three solutions of Gl, and each solution 
places a queen in row 0 at a different column. СТ interacts with the three relevant G® 
gadgets through the position of the queen on row 0. For example, if the queen is placed at 
(coll, 0) then the G® gadget for variable v; is forced to take its red solution. СТ does not 
directly interact with any other gadget. 

Definition 28 (Gadgets С? and G3). Suppose с = {u;,v;, к} where i <j < k, so 6; < 
9; < Ok. 


® (2 = (со + 1, С>, Ra, {}, {}) where: Cy = {со , coli}, and Ra = {0,7}. 
e (3 EI (colt + 1, Сз, Rs, {},{}) where: C3 = {colf, сой,}, and Rz = {0,14}. 


Gadgets С? and СЗ are illustrated in Figure 8. Both these gadgets have two solutions. 
They both interact with G° gadgets via the queen placed on row 0. They each force one 
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of the Сб gadgets to take its blue solution. С? and СЗ cannot both use the specific pair 
(coli, 0) because both have the same sum-diagonal. Taken together, С? and СЗ have three 
solutions, each of which forces two G° gadgets to take their blue solutions. 

It is natural to ask whether we could encode 3-SAT clauses using only Gt (adapted to 
use blue diagonals for the positive literals and red diagonals for negative literals), making 
С? and СЗ redundant. 3-SAT with at most 3 occurrences of each variable is NP-complete. 
However, a 3-SAT clause c is allowed to be satisfied by more than one of its literals. In 
this case the corresponding G! gadget would have more than one solution, so the reduction 
would not be parsimonious. Supposing we wished to prove only NP-Completeness we could 
use the simpler reduction from 3-SAT. However, we know of no way to simplify the other 
three reductions in Definitions 10, 16, and 20. 

The sum-diagonals for row 0 in Gl are д; + dee, 9; + cae Ok + Th os in each case a blue 
sum-diagonal from the relevant G° gadget. Similarly, the sum-diagonals for row 0 in G2, СЗ 
are д; + 4; ., д; + dic On + dies each a red sum-diagonal from the relevant G®°. 

The separations of 7 between rows and 21 between columns were chosen to make the 
following simple arithmetical fact true: 


Observation 29. Suppose that we have integers ал, а2, ст, со, 41, 42 with 0 < aı,aa < 2, 
0 < 41, 42 < 6, and suppose that 21(cı — cg) + 7(aı — a2) + (4: — 42) = 0. Then cı = cn, 
a, = а2 and dı = da. 


Proof. If cı # cg then |21(cı — с2)| > 21, while | + 7(aı — aa) + (а: — d2)| < (2- 0)7+ 6 = 
20 < 21. This makes it impossible to make the sum equal to zero unless cı = ca, and also 
7(aı — a2) (4: — 42) = 0. By similar reasoning, we must have aı = aa since if a; # aa then 
|7aı — 7ag| > 7, and |dı — da| < 6 < 7. Since we have cı = ca, aı = a2, we must also have 
dı = 42. 





























Before presenting the reduction itself, we prove (for all three gadget types) that the 
possible locations of the queens do not attack each other diagonally. 


Lemma 30. Within each of Gl, G2, СЗ, where the gadget has columns С; and rows Ву, the 
difference-diagonals a — В with a Е Ci, ВЕ В; are all distinct. The same applies for the 
sum-diagonals a+ В. 





Proof. We can deal with both cases at once by writing + where necessary. For v1,v2 € 
{i,j,k}, x € {b,r}, a1,a2 Е {0,1,2} we have a = coly, , В = 7aı. This makes a + В = 
dv, + di, с + 7aı. If two sum-/difference-diagonals within one G!, G?, or G are the same 
then we have (since values of x and с are fixed within a clause): 0 = 21(%1 — va) + 7(aı — 
az) + di, .— 4, с. Observation 29 gives 91 = v2 and aı = a2. In each gadget GL, G?, G the 
value of v and a uniquely determines a and В. Thus the value of a+ В uniquely determines 


a and В, which can only happen if all sum-/difference-diagonals are distinct. 


























Definition 31 (Reduction from Restricted 1-in-3-SAT to Problem 5). Consider an instance 
Mg = (V,C) of Restricted 1-in-3-SAT. For a clause с = {v;, vj, ик} € С we will write Ms, 
for the instance of Problem 5, Ms. = (GP; GF Gps Go, Ges Gahs For the reduction of the 
entire instance Mg we set: 


d 
Ms = LU M5, 
сЕС 
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Each G® can occur in the union up to three times, if v is in three clauses. However, 
each copy is identical since G° does not depend on с, so in the final set Mg, each G® occurs 
exactly once. 


Lemma 32. The possible sum-diagonals т Ms are all of the form dy +d,. + Та with 
0<a< 2, and additionally: а = 0 т G2; x =b in Gl; апах =r т G2, G3. Furthermore, 
if Ov, + але. + Taı = dy, + 4,2 с, + Таз, then we have all of v1 = v2, с1 = C2, 2 = 22, 
aı = a2. 


Proof. The first statement follows from a straightforward case analysis. In GO the possible 
values are exactly д, + dö.. In each of Gl, G?, G3 we have C C {col?, со, со; }, В © 
{7a | 0 <а < 2}, and so the sum-diagonal is col? + 7a = д, + &, + 7a. In Gj, « = b 
and in G2, СЗ, x = г, as required by the statement. For the second statement, suppose 
21(v1 — v2) + 7(aı — a2) + (dj, — 42 „,) = 0. Observation 29 applies to immediately give 

By construction of С, for a given variable 


U1 = V2, а1 = а2, and also dy!.., = doze. = Ay? ey: 
v1, the values с and x are uniquely determined from the value of dj, .. Hence cı = cg and 











v1 = %2. 





Lemma 33. For a clause с = {и оу, ок} from Me, Ms,c has exactly three solutions. In one, 
G? takes its red solution and GH, @ take their blue solutions; in a second G is red and 
G, a blue; and in the third a is red and G®, (© blue. 


Proof. We work by case analysis on where the queen in row 0 in G? is placed. This can be 
one of the pairs (col®, 0), (coll, 0), or (col®,0). In each case we will construct a solution set 
of queens S4(G) for each GE Мь,с, and show that 55. = {54(G) | СЕ Ms,-} is the unique 
solution to Ms. for that case. In the case for (со? , 0), the corresponding gadget G will use 
its red solution and the other two G® their blue solutions. 

First consider the case (col?,0) € S4(G!). This has the sum-diagonal col? = de. + 0j, 
the same value as a blue sum-diagonal of С9, so we set S4(G®) to its red solution. The 
other two pairs in $4(G!) must use the columns col! and col?. However, the pair (col?, 7) 
is excluded by D>. Hence we set $4(G!) = {(col?,0), (col!,7), (col? 14)}. Now consider 
G?. The columns available in row 0 are col? and col. But the pair (col, 0) has the sum- 
diagonal col; = dj, + 6;, and this is the same value as one red sum-diagonal of G°, and 
we are using the red solution for G® so we cannot use the pair (col?,0) in G2. So we set 
54(G?) = 1 (col},0), (col},7)}. But the sum-diagonal value of (col;,0) is the same as a red 
diagonal in GY, so we set S4(G¥) to its blue solution. Now consider СЗ. The pair (со! 0) 
is disallowed because it is on the same sum-diagonal as the pair (col7,0) € S4(G2). Hence 
we set $4(G?) = {(colf,,0), (coli, 14)}. The sum-diagonal value of (сой,, 0) forces us to set 
S4(G°) to its blue solution. 

The above has defined 55 completely for this case, and shown that there were по 
alternative choices at any point. This establishes uniqueness, but it is still necessary to 
show that 55,с is indeed a solution to Мс, i.e. that each S4(G) is a solution to G and that 
there are no sum-diagonal clashes between gadgets. The values $4(G®) are already known 
to be solutions. Each $4(G"), $4(G?), S4(G?) places exactly one queen in each allowed row 
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and column. Lemma 30 shows that no two queens can be placed on the same difference- or 
sum-diagonal within Gl, G?, СЗ. Also $4(G!) correctly does not contain the pairs (col?, 7) 
and (col? ‚ 14) which are disallowed by От. We now consider sum-diagonals from different 
gadgets. From Lemma 32, if two sum-diagonals have the same value of 6, + 4, ‚ + 7a, they 
have the same values of v,c,x,a. We only have to consider а = 0, since from Lemma 32, 
СУ only has a = 0, С? and СЗ do not share a non-zero value of a, and G! does not share a 
value of x with G?, СЗ. For a = 0 the only cases of duplicate values of д, + 4, с are: between 
СО and Gt, С2, G?; and the value 6; + dj, between С? and G3. The construction of 55 с 
above avoided using any of these sum-diagonals twice, and thus is indeed a solution. 

We omit details of the cases (col, 0) and (col?,0) in Gt, which force Gy and G® (re- 
spectively) to take their red solutions with the other copies of С° blue. The case (col? 0) 
is exactly similar to the case above except that Оу rules out the pair (col?®, 14). There is a 
slight difference for the pair (coll, 0) in G!: the solution of G! is fixed by the exclusion of 
both diagonals in От, while the solutions of G? and СЗ are unique because (coli, 0) cannot 














be used in either since it has the same sum-diagonal as the red solution to GH. 


Lemma 34. For any truth assignment Sg which solves Mg then 55 is a solution to Ms, 
where 55 = {54(G) | G € Ms}; and where Sı(G®) is the red solution iff S6(vi) = true and 
the blue solution otherwise; and where values of S4(G'), Sı(G?), Sı(G?) are set as in the 
proof of Lemma 33. 


Proof. Because Sg is a solution to Mg = (У, С), for any clause с = {v;,v;,up} € С we know 
that Sg assigns exactly one of the literals in с to true. Hence exactly one of 54(С°) is red 
and the other two are blue, and, by Lemma 33, the values of $4(G!), 54(С2), Sı(G?) are 
well-defined solutions and unique. It remains to show that for any two gadgets G,G’ € Ms, 
and any two queens а # а’, q € S4(G), q € S4(G’), the sum-diagonals of а and 4’ are 
different. Lemma 32 shows that if q and а’ have the same sum-diagonals, then they have the 
same values (in terms of that lemma) of v; and va, and of cı and со, that is they arise from 
the same combination of clause and variable. Hence either G = С" = G® ог С, С" € M;.., or 
both. In either case, their sum-diagonals must be different according to Lemma 33. Finally 
we note that $5 is well defined and unique because each 54(С°) is forced to be either the 
red or blue solution by the truth value of v;, even if v; appears in no clauses. 














Lemma 35. For any solution 55 = {S4(G) |G € Ms} to Ms, then 56 is a solution to Мб, 
where 56 (%) = true iff S4(G®) is the red solution and Sg(v;) = false iff Sı(G®) is the blue 


solution. 


Proof. Sg is well defined since each G? can only take its red or blue solution. From Lemma 33, 
exactly one variable from each clause in Mg is set to true by Se. 














Theorem 36. The reduction of Definition 31 from Restricted 1-in-8-SAT to Problem 5 is 
a parsimonious reduction. 


Proof. From Lemma 34 and Lemma 35, and the fact that the maps defined in those lemmas 
are inverses of each other. 














Theorem 37. The n-Queens Completion problem (Problem 2) is NP-Complete and #P- 
Complete. 
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Proof. Theorems 36, 23, 19, and 13 give successive parsimonious reductions starting from 
Restricted 1-in-3-SAT and ending with n-Queens Completion. By inspection each reduction 
can be computed in polynomial time. Restricted 1-in-3-SAT is NP-Complete and #P- 
Complete by Theorem 4 so this gives NP- and #P-hardness. NP- and #P-easiness follow 
since potential solutions to n-Queens Completion are straightforward to check. 














As an intermediate problem, the Excluded Diagonals Problem (Problem 4) is also NP- 
and #P-Complete. As a straightforward corollary we obtain the NP-Completeness of the 
“Blocked n-Queens” problem, which has been used as a benchmark for answer set program- 
ming (ASP). In this problem, some squares are blocked from holding queens (Namasivayam 
& Truszczynski, 2009). 


Problem 7 (Blocked n-Queens). PROBLEM: М7 = (n, P) where т is an integer and P is 
a set of queens (possibly attacking each other) that all fit on board size п. SOLUTION: А set 
57 of queens which is a solution to the n-Queens Problem рот n and such that PN S7 = {}. 


Corollary 38. The Blocked n-Queens problem is NP-Complete and #P-Complete. 


Proof. The reduction from n-Queens Completion to Blocked n-Queens is straightforward. 
For a problem Ma = (п, Pa) we block each queen on the same row as one of the queens in 
P,. That is, we form the set Р; = Ua Bert B)|O<i<n,i#a}. Any solution to M7 
= (п, Pz) necessarily contains each queen in Pa so is a solution to Ma. This establishes NP- 
and #P-hardness, while -easiness is true because potential solutions to Blocked n-Queens 











problems can be checked easily. 





5. Empirical Study 


We have shown that n-Queens Completion and Blocked n-Queens are NP-Complete, but 
how hard are they to solve in practice? This is an important question given the repeated 
controversy surrounding the use of n-Queens as a benchmark problem discussed in Sec- 
tion 1. We first describe solvers based on three different technologies: a special purpose 
solver; a SAT solver; and a constraint solver. We describe methods of generating random 
benchmark instances of three problem classes: Blocked n-Queens, n-Queens Completion, 
and the Excluded Diagonals Problem (Problem 4). We looked for phase transitions (Hogg, 
Huberman, & Williams, 1996) between solvable and unsolvable instances, in common with 
other NP-Complete problems, for example Boolean satisfiability (SAT) (Mitchell, Selman, & 
Levesque, 1992; Achlioptas, 2009; Altarelli, Monasson, Semerjian, & Zamponi, 2009). It is 
common in NP-Complete problems to see that randomly generated instances are hardest at 
the transition from solvable to unsolvable instances (Cheeseman, Kanefsky, & Taylor, 1991; 
Mitchell et al., 1992; Gent, MacIntyre, Prosser, & Walsh, 1995). We do see this pattern 
for every solver in two of our generators but not in the third: a straightforward generation 
method for n-Queens Completion appears unable to produce instances that are hard to solve 
using SAT. 

Experiments were performed on a 32-core AMD Opteron 6272 at 2.1 GHz with 256 GB 
RAM. Savile Row was executed in OpenJDK Java 1.7.0_ 131 with the 64-bit server VM. For 
random instances, we solved 1000 randomly generated instances per data point, and solvers 
were allowed to run to completion. We have made available our code including the adapted 
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version of Savile Row that we used,’ detailed results including some additional graphs,® and 
experimental instances.? 


5.1 Three Solvers 


We constructed three types of solver, each of which is complete, i.e. is guaranteed to either 
find a solution or confirm that there are none. We make no claim as to how good these 
methods are compared to other possible methods: we wrote solvers using three different 
technologies to ensure that if we found hard instances, they were hard for at least a variety 
of solvers. 

The first approach is a special-purpose solver written in C++. It is based on bit manip- 
ulation techniques previously used for counting solutions of n-Queens (Richards, 1997). It 
searches exceptionally fast per search state although the intelligence of the search process 
is limited. The algorithm is a backtrack search that constructs a solution starting at row 
0 and moving up the board. It uses three bit patterns to represent attacks on the current 
row В via the sum-diagonals (148), columns (colsg) and difference-diagonals (rdg). Each 
bit pattern contains a 1 in position a iff a queen attacks position (a, 3) via the respective 
diagonal or column. When a queen (a, В) is placed on the board, bit « is set in Idg, colsg, 
and rdg. When moving from row В to 6+1, ldg41 = в << 1 (i.e. Id is shifted one position 
left), rdg41 = rdg >> 1 (rd is shifted one position right) and cols3,1 = colsg. The attacked 
positions of row В + 1 are given by the bitwise OR of в-1, colsg41, and rdg41. Another bit 
pattern technique is used to find the rightmost non-attacked position that has not yet been 
tried. We made simple adaptations to account for the fact that some queens are already 
placed in n-Queens Completion. First, we skip rows that contain a queen. Second, ldo, 
colsy, and rdg are populated with attacks of the pre-placed queens. Third, when shifting 
ld (rd) the least (most) significant bit is populated from an array to allow for a pre-placed 
queen attacking diagonally from above. A different adaptation was made for Blocked n- 
Queens where the blocked positions of a row В are masked out using a bitwise AND before 
attempting to place a queen on row В. А third adaptation was made to solve the Excluded 
Diagonals Problem, in which [4% and rdp are pre-populated with excluded diagonals, and 
when shifting Id (rd) the least (most) significant bit is populated from an array containing 
the excluded diagonals that do not intersect row 0. 

The second approach is to encode to SAT and apply a CDCL SAT solver. We use the 
SAT solver Lingeling (Jarvisalo, Heule, & Biere, 2012) which was winner of the Sequential, 
Application SAT+UNSAT track of the SAT competition 2014.19 We wrote two Python pro- 
grams to encode n-Queens Completion and Blocked n-Queens into SAT. SAT variables are 
added for each square which is not attacked by a preplaced queen (for n-Queens Completion) 
or not blocked (for Blocked n-Queens). We use the commander-variable encoding (Klieber 
& Kwon, 2007) for Boolean sum < 1 constraints (using a group size of 3) to add ‘exactly-one’ 
constraints for each row and column, and ‘at-most-one’ constraints for each diagonal. To 
solve the Excluded Diagonals Problem we implemented the reduction to n-Queens Comple- 
tion (Definition 10 and Definition 16) with a separate Python program. While this approach 





7. https://ipg.host.cs.st-andrews.ac.uk/n-queens-completion-experiments-code.tgz 

8. https://ipg.host.cs.st-andrews.ac.uk/n-queens-completion-experiments-results.tgz 
9. https://ipg.host.cs.st-andrews.ac.uk/n-queens-completion-experiments- instances. tgz 
10. 


Lingeling version ayv 86bf266b9332599f1b876e28a02fe8427aeaa2db 
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Figure 9: Random Blocked n-Queens instances: y axis shows probability of a solution ex- 
isting (left) and mean nodes searched by Lingeling (right) against x„,, on x axis. 
The colour and plot symbol indicate the value of n. 


generates large n-Queens Completion instances, the number of non-attacked squares (and 
therefore the SAT encoding) is small. Further, since we can generate hard instances of the 
Excluded Diagonals Problem, this approach gives us confidence that we can also generate 
hard instances of n-Queens Completion. 

The third approach uses the constraint solver Minion 1.8 (Gent et al., 2006), with a 
constraint model optimised by Savile Row 1.6.4 (Nightingale et al., 2017). For each row 
В, the model has three variables representing the column, sum- and difference-diagonal 
of the queen placed on row В. The model exploits powerful constraint propagation by 
using three all-different constraints to enforce non-attacks on columns, sum- and difference- 
diagonals. Two table constraints are used on each row to connect the column variable to the 
other two variables. Savile Row preprocesses each instance with Minion which may assign 
variables. Then it applies partial evaluation to simplify the constraints before generating the 
Minion input file. n-Queens Completion is implemented by assigning the appropriate column 
variables, Blocked n-Queens by removing values of the column variables, and Excluded 
Diagonals Problem by removing values of the variables representing diagonals. 


5.2 Blocked n-Queens 


We start with Blocked n-Queens (Problem 7) since it has previously been used as a bench- 
mark for ASP solvers (Namasivayam & Truszczynski, 2009; Gebser, Liu, Namasivayam, 
Neumann, Schaub, & Truszczyriski, 2007). We proved the problem NP-Complete and #P- 
Complete in Corollary 38. The competition instances came from a set of 100 for each 
n € {28, 48, 50,56}.!! The SAT solver was able to solve all instances within 8 minutes CPU 
time, the hardest taking 431.4s total time and 2,890,435 nodes. 

To investigate further we designed our own generator, since the generation method for 
the competition instances is not published. A row or column with zero unblocked squares 





11. Thanks to Martin Gebser for providing the 400 instances. 
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makes an instance trivially unsolvable. Nine of the ASP competition instances (5 at п = 56) 
are trivial for this reason. More subtle problems remain with up to 3 unblocked squares. 
For example, if the only three unblocked squares in row 1 were in columns 5, 6, and 7, 
then row 2 could not have a queen in column 6. This is an implicit unary constraint 
which might combine with others to make a problem trivial: similar issues with ‘flawed’ 
instances have been seen in generation methods for several other combinatorial problems 
(Achlioptas, Molloy, Kirousis, Stamatiou, Kranakis, & Krizanc, 2001; Gent & Walsh, 1999; 
Gent, Nightingale, Rowley, & Stergiou, 2008). To avoid this danger, we require each row 
and column to have at least 4 free squares. We generate an instance with b blocked squares 
as follows. We start with an empty set of blocked squares. We select a non-blocked square 
uniformly at random, rejecting any that leave a row or column with 3 or fewer unblocked 
squares, and add it to the set. Squares are added to the set until we have 6 blocked squares, 
or all remaining non-blocked squares are rejected (in which case we restart with an empty 
set). 

We calculated а constrainedness parameter к (Gent, MacIntyre, Prosser, & Walsh, 1996). 
The expectation is that there should be a phase transition at a critical value of « < 1, with 
search hardness peaking near that value. For a random ensemble with state space 5 and 
probability р; that a given state se 5 is a solution, к = —1055 р. /1055 |S|. We take 5 
to be the set of solutions to the n-Queens problem. It is conjectured that |5| ~ n!/2.54” 
(Cloitre, 2002): by Stirling’s formula, |5| = V2rn(n")/(2.54e)"”. With b blocked squares, 
each queen in s € 5 has probability b/n? of being disallowed, so ps = (1 — b/n?)”. Putting 
the approximations together and simplifying, we get: 





a, п(оват — log,(n - b/n)) is 
> n(logy n — logy 2.54e) + $(1+ logy п + logy п) 


Figure 9 shows experimental results using the SAT solver. We see a transition in solvabil- 
ity at a critical value of к = 0.97, with the transition sharpening as n increases. There is a 
hardness peak associated with the solvability transition in each solver. We only show results 
for the SAT solver as it is much better than the other two, being able to solve instances 
up to п = 60 in reasonable time (peak mean time 681.9s, peak mean nodes = 3.5 x 10°). 
The Minion-based solver solves only up to 40 (peak mean time 358.9s, peak mean nodes 
= 1.8 x 10°), and the special purpose solver up to 35 (peak mean time 126.7s, peak mean 
nodes = 2.6 x 10°). Increasing п, by 5 above these values made the experiments prohibitively 
expensive to run with expected peak mean CPU times above 1000s in each case. Examining 
Equation (1), the only influence that b has is in the term п — b/n, which is a natural param- 
eter, i.e. the mean number of unblocked squares in each row. Given the observed critical 
value к = 0.97, solving for п — b/n shows that аз п increases the critical value of п — b/n 
grows very slowly, ~ (2.54е)9.97п0-03 ~ 6.50.93. Figure 10 shows our experimental results 
for SAT plotted against п — b/n instead of ky. Results are consistent with this analysis 
but are over too limited a range to comment on the scaling. 


5.3 Random Placement of Queens 


The most natural way to generate random instances of n-Queens Completion is to place 
queens at random positions. Given n and a number p of pre-placed queens, we start with 
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Figure 10: Random Blocked n-Queens instances: у axis gives probability of a solution ex- 
isting (left) and mean nodes searched by Lingeling (right) against n — b/n (mean 
number of unblocked squares per row) on x axis. The colour and plot symbol 
indicate the value of п. 


an empty list P = {} of pre-placed queens, adding queens iteratively until |P| =p. At each 
step a queen q is chosen uniformly at random from the set of queens that do not attack 
any queen already in P. If we allowed queens in P to attack each other then the birthday 
paradox would mean that we could only generate instances with р = O(,/n) before almost 
all instances contained attacks: this would be a flawed generator. We typically experimented 
with fixing the difference п — р, the number of queens to be placed. We do not show detailed 
results here but graphs and detailed results are included in our online archive. Results 
with the special purpose solver were encouraging. We saw a transition between solvable 
and unsolvable regions, associated with hard instances. The hardest instances can require 
an astronomical number of search nodes: one instance at n — p = 30,n = 108 required 
= 6.64 x 10!? nodes, search taking = 1.8 days despite = 42 million nodes being searched 
per second in a single thread. There was also an apparent hard region with very small 
values of p, similar to clusters of difficult instances in the mainly satisfiable region in other 
problems (Hogg & Williams, 1994; Gent & Walsh, 1994). However, when using the SAT 
solver there are very few difficult instances. Brief investigation of much larger instances (e.g. 
n — p = 500) failed to find any the SAT solver could not solve in a short time and even 
instances with n in the millions had solutions. While not definitive, these results suggest 
that this generation method does not lead to instances that are intrinsically hard, and it 
may not have a phase transition at all for large n. Therefore we conclude it is not a suitable 
generator for n-Queens Completion benchmarks. 


5.4 The Excluded Diagonals Problem 


While the simple and natural generator for random n-Queens Completion seems to be unable 
to generate hard instances, we found that the following generation method for Problem 4 
(the Excluded Diagonals Problem) generates consistently hard instances for the three solver 
types we consider. Given n, we allow all rows and columns but select d different diagonals to 


841 


GENT, JEFFERSON, & NIGHTINGALE 






































© © 
- | о 
ЗЕ Ен: 
Ф a EB 
Er, 
— ® 
5 © at Е 
$ | PLUS ISCH Re 
- BER o> BER LIE: 
Oe eu Е 8 og PU Fs 
5 + 988° ТО ОЗУ уу" ое 
$ 2.93% SEES See SO el “$e ga 
+ 2 Е ЗАРЕ $ 
5 1 АД 
< io 9-0. AAO 
м | ö 
© - 
о 3 
> | rape 
= 2 т T т 
0.0 0.5 1.0 1.5 
© ioe] 
© © 
+ 4 7 4 
Ф Ф 
© © Bea 
о © ве Ва 
57 6 1 FLW Е: 
u 2 20,0 ФФ оби m 
: | cqugeatttensseescea gett 
3 5 BEERS SO ФФ? oy 
ci т ak & So BB $ on 
2 | geete re ов 
RG XXX XX X OSV GS $ 
a a 3 $ + Xx Cove op 
= ¢ 4 O0 AT Ex Love adh 
© © ie = 0-0.A KR FS 
$ $ о 
ATI SEN 
> © Os KAKO SS ee 
aa eal eB 
oy | | | = u | | —— 
0.0 0.5 1.0 1.5 0.0 0.5 1.0 15 


Figure 11: Random instances of the Excluded Diagonals Problem. The x axis represents d/n 
on all plots, and the colour and plot symbol indicate the value of n. Probability 
of a solution existing on y axis (upper left), mean nodes searched using special 
purpose solver (upper right), Lingeling (lower left) and Minion (lower right). 
The key is the same for all plots and the range of the у axis is the same for all 
hardness plots. 


exclude out of the 4n—2 possible diagonals (i.e. 2n—1 each of sum- and difference-diagonals). 
That is, we construct an instance Ма = (n,{},{},D~,D*) subject to |D~|+|Dt| = а. 
Diagonals are selected one by one, randomly from those remaining, but rejected if excluding 
the diagonal leaves any row or column with 3 or fewer unblocked squares. For d > 2n — 2 
instances are trivially unsolvable: a solution requires 2n allowed diagonals. 

Results with our three solver types are shown in Figure 11. We do observe a transition in 
solvability. It appears to be sharper as n increases. For these sizes the transition in solvability 
occurs at d/n = 1, though we cannot say if that trend continues as n increases. All of the 
solvers show a hardness peak associated with the solvability transition. It is interesting that 
solving is so hard that we could only experiment up to n = 21. The maximum mean time 
taken for the constraint solver was 262.7s, 78.3s for the SAT solver, but only 4.6s for the 
special purpose solver. The special purpose solver searches many times more nodes than the 
more intelligent methods, but it searches so fast that for these instances it is the fastest. 
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The Excluded Diagonals Problem may be translated to both n-Queens Completion and 
Blocked n-Queens so can be used to provide hard random instances of both problems. Some 
interesting questions remain for this generator. We do not know if this generator continues to 
produce hard instances for larger n, and if so where the phase transition point is. Although 
the instances can easily be translated to Blocked n-Queens, the definition of к above is 
not applicable here because the resulting instances are not drawn from the same random 
distribution: it would be interesting to find an appropriate definition of «. 


6. Conclusions 


We have shown that the n-Queens Completion problem, the Blocked n-Queens problem, 
and the Excluded Diagonals Problem are NP-Complete and #P-Complete. For anybody 
who understands the rules of chess, n-Queens Completion may be one of the most natural 
NP-Complete problems of all: “Given an n x n chessboard on which some queens are already 
placed, can you place a queen in every remaining row so that no two queens attack each 
other?” 

We have presented generators for hard random instances of the Blocked n-Queens and 
Excluded Diagonals problems, shown that this hardness persists across three solvers using 
very different technologies, and that the hardness is associated with a phase transition in 
solvability. A natural model for n-Queens Completion does not appear to generate consis- 
tently hard instances. We have also shown that care must be taken in generation methods 
to avoid flaws which can be found in extant benchmarks. Our experiments raise many inter- 
esting questions for future research. These include the asymptotic scaling of the transition 
in Blocked n-Queens and the Excluded Diagonals problems. An interesting open problem is 
to find a direct generator of hard random instances for n-Queens Completion which works 
by randomly placing queens on a chessboard without using reductions. 

The importance of our work is that the n-Queens problem has been a very widely used, 
but often criticised, benchmark in Artificial Intelligence. Our work shows that the n-Queens 
Completion problem or the closely related Blocked n-Queens problem can serve as valid 
benchmark problems. We implemented three complete methods using different underly- 
ing technologies. The best of these in our experiments is to translate into SAT using the 
commander encoding and solve with a modern SAT solver, but it would be interesting to 
compare this with more of the very many techniques that have been applied to the n-Queens 
problem. Solvers for n-Queens typically require minimal adaptation to work on the variants 
in this paper. By providing random instance generators we hope to encourage researchers to 
test their n-Queens algorithms on variants that can better discriminate between good and 
bad solvers than the original problem. 
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